Ab initio study of ladder-type polymers polythiophene and polypyrrole 
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This article presents an ab initio study of four polymers, polythiophene, polypyrrole, ladder-type 
polythiophene, and ladder-type polypyrrole. Upon an analysis of the variation of the band gap 
when comparing the unconstrained and the ladder-type polymers, a discrepancy was found between 
the thiophene and the pyrrole polymer families. For polythiophene, the ladder-type polymer has a 
larger gap than the unconstrained polymer whereas the opposite is found for the pyrrole polymers. 
The structural properties and the charge densities using the Bader charge analysis of these four 
compounds are investigated. The different band gap behaviors in thiophene and pyrrole polymers 
can be explained in terms of the competition between the bond length alternation and the effect of 
the charge density in the carbon backbone. 



INTRODUCTION 

Conjugated polymers raise much interest because of 
their intriguing electronic properties and their foreseen 
technological applications. They have many advantages 
when compared to inorganic semiconductors such as their 
easy processing, their tunable optical gap and their fa- 
vorable structural properties. Their electronic properties 
are mainly due to their delocalized 7r-electrons along their 
carbon backbone. Polyanniline 1] 2], polypyrrole [3j and 
polythiophene are examples of polymers with promis- 
ing future developments. In particular, polythiophene 
and its derivatives were used in several applications 
such as numeric display device Q, surface light emitting 
diode(SLED) and light emitting diode(LED)@. 

One particular class of polymer that is of present inter- 
est is the ladder-type polymers. These have additional 
bonds as compared to other polymers which link the 
neighboring monomers rigidly together eliminating the 
possible dihedral degree of freedom. These polymers are 
known generally to exhibit small band gaps, due partly 
to their planar configurations which maximize the align- 
ment of the 7r orbitals. Moreover, these ladder-type poly- 
mers have the potential to exhibit very high intrachain 
mobility 0- 

Recently, Oyaizu and al. || reported the first synthe- 
sis of a ladder-type polythiophcnc(LPT) and published 
details about its electronic structure that they have char- 
acterized both theoretically and experimentally. Using a 
parameterized Hartree-Fock model (PM5), they found a 
band gap reduction in LPT as compared to polythio- 
phene(PT); they also obtained indirect experimental ev- 
idences of this gap reduction. They relate this behavior 
to a reduction of the bond length alternation observed in 
their calculations when the backbone of the LPT is com- 
pared to PT. This argument is similar to the one used 
in polyacetylene which explains the band gap and the 
dimcrization of the atomic structure. 

In this article, a pseudopotential density-functional 



theory study of the electronic properties of four polymers, 
the polythiophcnc(PT), ladder polythiophene (LPT), 
polypyrrole (PPy), ladder polypyrrole (LPPy) is pre- 
sented. The results indicate a decrease of the band gap 
for the ladder- type version of polypyrrole, however, on 
the contrary, the ladder-type polythiophene present a 
larger band gap as compared to PT. This difference be- 
tween these two systems is examined in detail in the re- 
maining of this article. 



COMPUTATIONAL METHODS 

The results reported in this article were computed 
within the framework of the density-functional the- 
ory(DFT) as implemented in two codes that use differ- 
ent basis set to represent the electronic degree of free- 
dom. The first is the Abinit package^ which repre- 
sents the electronic states with a plane- wave basis set, 
useful in the pseudopotential formulation of periodic 
structures. The exchange-correlation energy was cal- 
culated both in the local density approximation(LDA) 
using the Teter-Pade parametrization which reproduces 
the Ceperley-Alder data|Kj and in the Perdew-Burke- 
Ernzerhof(PBE) generalized-gradient approximation llj. 
The pseudopotentials were generated with the Trouiller- 
Martins scheme^^ an d their portability was fully tested. 
Numerical convergence of the total energy within 1 
mHa/atom was reached for a number of plane- waves in 
the basis set corresponding to a kinetic energy cutoff of 
35 Ha and a sampling of the Brillouin zone of 8 k-points 
on a shifted grid. 

The Abinit package is a solid-state oriented code, 
which implicitly generates Born-von Karman periodic 
boundary conditions in all directions. This implies that 
within Abinit, the electronic structure of a single polymer 
is not directly computed, but rather the electronic struc- 
ture of a uniform array of polymers is simulated. The un- 
wanted interaction between neighbor polymers can easily 
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be minimized by increasing the space between polymers 
which however, also augments the computational efforts. 
This approach has the advantage that the electronic de- 
grees of freedom are represented by an orthogonal basis 
set which covers all space equivalently and that can be 
check for completeness simply by increasing the number 
of plane-waves, e.i. the kinetic energy cutoff mentioned 
above. 

The charge densities presented in Fig. ^ was obtained 
using Bader charge density analysis 13] . To get con- 
verged results, a basis that includes plane-waves up to a 
kinetic energy of 100 Ha was required to get the charge 
density on a very fine grid in real space. 

Furthermore, to compare the validity of the function- 
als used and of the results for isolated polymers, the 
Gaussian 03 codeQ was employed. This code is ori- 
ented towards molecules and represents the wave func- 
tions with a gaussian basis set. The basis used was 6- 
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will be explained later. Within Gaussian 03, periodic- 
ity can be imposed only in the directions of interest. In 
the present case, only one translation vector is specified 
to obtain truly one isolated polymer. The exchange and 
correlation functionals used were again the LDA and the 
PBE to compare the results for a single polymer. As 
can be seen in Table the band gap values obtained 
with Abinit and Gaussian for the PPy and LPPy are 
similar, but for the PT and LPT a small difference is 
noted. This variation can be explained, in part, by the 
different pseudopotentials used in the two simulations. 
Nevertheless, the band gap variations reported in Ta- 
ble [H] for these polymers are comparable and differ only 
by 0.1 eV between the Abinit and the Gaussian results. 
To investigate the validity of the bond length alternation 
obtained with LDA and PBE, the calculations with the 
B3LYP functional[24| was also employed which includes 
exact exchange. This functional is implemented in the 
Gaussian 03 code but not in the Abinit package. Finally, 
the Gaussian 03 code used an equivalent of 96 k-points 
to integrate over the Brillouin zone. 

The polymers have been initialized in a coplanar con- 
formation and have been fully optimized. All the poly- 
mers are considered uncharged. Although the relaxation 
of the atomic structure was not constrained the coplanar 
symmetry was never broken during the optimization pro- 
cess, both with the Abinit and the Gaussian codes. This 
is consistent with the expectation of an uncharged ir- 
conjugated system and with the results obtained for the 
oligomer versions of these polymers [2^. The treatment 
of periodic systems within the Gaussian code is a fairly 
new functionality. Unfortunately, the B3LYP functional 
is not well suited for periodic systems. The long-range 
part of the exact exchange in this functional causes con- 
vergence problems which are correlated to the basis set. 
It was found that the 6-31G basis set was the largest 
that will converge for all the polymers. This correlation 
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FIG. 1: Atomic structure of the polymers and oligomers. The 
charges and bond length obtained within LDA are indicated 
on the atomic structures of the polymers. 

is simple to understand, reducing the basis set effectively 
reduces the long-range overlaps between wave functions 
making the exact exchange part of the B3LYP functional 
more manageable. Fortunately, these long-range effects 
have only little influences on the properties studied and 
the 6-3 1G basis set is sufheiant the treat the present elec- 
tronic systems. 

RESULTS AND DISCUSSION 
Thiophene based polymers 

The first synthesis of LPT was reported by Oyaizu and 
al.Q. They explored the electronic structure of this poly- 
mer with indirect experiments and by performing calcula- 
tions using the PM5 level of theory, which uses empirical 
data to parametrize the two center integrals within the 
Hartree-Fock theory. This method explicitly includes ex- 
change contribution, but partially includes correlation ef- 
fects through the empirical fitting of the electronic prop- 
erties. Within this method, the atomic structure of the 
PT is characterized by a clear bond length alternation 
in the carbon backbone. The bond opposite to the sulfur 
atom within the thiophene monomer is practically the 
same length as the intermonomer bond, being only 0.011 
A shorter, whereas the other C-C bonds within the thio- 
phene is 0.073 A shorter. The bond alternation obtained 
is therefore about 8r ~ 0.06 A for PT. However, with 
the LPT atomic structure, there are only two inequiva- 
lent C-C bonds and they are found to be roughly of same 
length resulting in a bond alternation of just Sr = 0.005 A 
for this polymer. All carbon atoms are equivalent in the 
LPT, whereas there are two inequivalent carbon atoms 
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FIG. 2: LDA band structure of PT and LPT 



in the PT. 

The authors then discuss the electronic structure ob- 
tained for the optimized structure. The band gap for 
both polymers corresponds to a direct transition at k=0. 
The calculated band gap for PT is found to be 6.45 eV 
while for LPT it is found to be 5.92 eV. They argue 
that the missing correlation effects should only result 
in a constant shift of the unoccupied electronic levels, 
which should reduce the band gap of PT to the exper- 
imental value of 2.1 eV. Applying the same shift leads 
to a band gap of 1.5 - 1.6 eV for LPT. The reduction 
of the band gap between the two polymers is then ex- 
plained using the reduced bond alternation by a tight- 
binding-likc argument as it is done in the polyacetylene 
case. This result is compared with experiment via UV 
absorbtion measurements, photolumincsccncc measure- 
ments and electrochemical properties (redox reactions) 
of thin films. The main difficulty of these procedures are 
that LPT is insoluble, forming a opaque film, and cannot 
be easily functionalized. This means that direct measure- 
ment of the UV absorption and the photoluminescence of 
LPT in solution is not possible. Results for the oligomers 
and the closest soluble precursor of both polymers are 
presented and suggest a band gap diminution. They also 
report results on both non-ladder oligomers and ladder 
oligomer. Fig. ^ shows the non-ladder trimer oligoth- 
iophene(DS-T 3 ) and its ladder version(Me 2 FT 3 ). The 
UV absorption spectra and the photoluminescence inten- 
sity maximums are both found to be red-shifted for the 
ladder version of the oligomer. 

The results obtained in the present study are in con- 
trast with the results presented above, the band gap of 
LPT is found to be 0.18 to 0.33 eV larger than PT de- 
pending on the functional used, as shown in Tables U 
and [H] The LDA band structure of these two poly- 
mers are depicted in Fig. |3 The band gap obtained 



with B3LYP is much more consistent with the experi- 
mental value of 2.1 eV for Pt[2(|. After optimization of 
the structures with the B3LYP (LDA) functional, the 
bond length alternation is found to be Sr ~ 0.04(0.03) A 
for PT and Sr = 0.016(0.003) A for LPT. The B3LYP 
results show that the bond alternation is increased com- 
pared to LDA and PBE functionals, which as already 
been seen in the case of trans-polyacetylene, thiophene 
e F based oligomers^ and polyyne oligomers [27|. 
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TABLE I: Calculated band gaps for different functionals of 
the selected polymers given in eV. 
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0.335 


0.309 
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0.179 


LPPy - PPy 


-0.182 


-0.24 


-0.080 


-0.078 


-0.231 



TABLE II: Variation of the band gaps of Table □ given in eV. 

The primary concern will be the variation of the 
band gap between a polymer and its ladder-type equiv- 
alent. The net result of the increased bond alternation in 
B3LYP is an almost uniform increase in the band gaps 
of all the polymers. Therefore, the difference between 
B3LYP and LDA functionals will not be an issue since 
both give the same variational behaviour, as illustrated 
in Table ITU for clarity. 

A survey of the literature reveals that the standard be- 
havior as been a decrease of the band gap in the ladder- 
type form of a polymer. It is the case for benchmark 
polymers like polyfluorene and polycarbazole[2^. The 
band gap reduction in these polymers can be readily un- 
derstood by the structural changes incurred during the 
ladder transformation. The unconstrained polymers are 
not coplanar: they have a non-zero value for the dihe- 
dral angles between monomers, varying between 26 and 
27 degrees which is a result of the electrostatic repulsion 
between the hydrogen atoms or other side groups. This 
causes a misalignment of the tt orbitals that decreases the 
dispersion of their bands resulting in an increase of the 
band gap. On the other hand, the planarity is enforced by 
the additional bonds in the ladder- type polymers. The 
7r orbitals interaction between monomers is maximized 
which favors a delocalisation of the wave functions and 
cause an increase in the dispersion of these levels, hence, 
a reduction of the band gap. 

Generally, the ladder-type polymers exhibit smaller 
band gaps then their unconstrained counterpart because 
they are planar. However, this cannot be the case 



for polythiophene, since the original polymer is alreai 
coplanar. The variation of the band gap must therefo 
be linked to more subtle changes of the geometry. 

As noted by Chung and al. 29], the absorption propc 
ties of PT closely resemble that of irans-polyacetyler 
The backbone consists of trans segments linked throuj 
a cis-like unit which implies that its electronic pro 
erties will lie between the two structures. A simp 
nearest-neighbour tight-binding model of the 7r-electro 
of <rans-polyacetylene shows that the bond alternatii 
will be critical in the determination of the band ga 
All the carbons atoms (with their corresponding hydr 
gen) in this structure are equivalent. Hence, if the bo: 
lengths are equal, the matrix elements corresponding 
tight-binding hopping parameters are also all equal. Tl 
leads to a closure of the band gap at the Brillouin zo: 
edge. As mentioned before, in the case of PT the ce 
bon atoms in the backbone are not all equivalent, whi 
suggest that even if the bond alternation is suppresses, 
the band gap should decrease but might not completely 
vanish. 

The key assumption in the preceding argument was 
the decrease in the bond alternation while preserving 
the ircms-polyacetylene like conformation. However, in 
the ladder- type polymers, the backbone has changed as 
compared to the non-ladder version and they now re- 
semble more closely to that of cis-polyacetylene. In the 
cis case, the carbons atoms are again all equivalent, but 
the hopping matrix elements are not. This is a subtle 
consequence of the geometry. A close consideration of 
the position of the hydrogen atoms reveals the change in 
their orientation between two consecutive pairs of carbon 
atoms. Thus, even with equal bond lengths there is no 
symmetry enforcing the values of the hopping terms to 
be equals and the band gap at the zone edge does not 
vanish. 

The present results can now be understood with these 
arguments. The ladder- type polymer does decrease the 
bond alternation in the backbone, but in the new con- 
formation, this does not reduce the gap. A perfect ex- 
ample of this is the 1.148 eV band gap with no bond 
alternation found in LPT with the LDA formalism ob- 
tained with the Gaussian code. The main reason for this 
band gap variation is due to the changes in the atomic 
structure between the non-ladder and the ladder types 
polymers. In the non-ladder type, the polymer is made 
of trans- and czs-segments whereas the ladder-type is en- 
tirely made of cis-segments. Care must be taken when 
comparisons are made with experiments, because the re- 
sults stated for the polymers will not necessarily be cor- 
respond to the oligomers. The tight-binding model and 
the band representation, breaks down without the Born- 
von Karman periodic boundary condition which means 
that a direct generalization to the oligomer cases cannot 
be made. Even the electronic structure of the precur- 
sors may be dominated by others effects related to the 
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FIG. 3: LDA band structure of PPy and LPPy. The dashed 
line represents the nearly free electron states(NFE). 

side-chains, as explained later in the case of PPy. 

To address the experimental results directly, the 
HOMO-LUMO separation of both oligomers showed in 
Fig. n using the B3LYP functional were calculated. A 
value of 3.25 eV is obtained for DS-T3 and 3.53 eV for 
Me2FT3. For simplicity, the side-chains in DS-T3 were 
not explicitly included since they are not expected to 
influence the HOMO or LUMO energies. Even in the 
case of the oligomers the band gap is found to increase 
in the ladder-type conformation which is contrary to the 
reported behavior observed in the optical spectra. The 
differences may be the consequence of a solvent inter- 
actions or of exciton effects that are not fully taken in 
consideration in the present calculations. 

Pyrrole based polymers 

Interestingly, the simple tight-binding model seems to 
fail for the polypyrrole family. All the DFT calcula- 
tions shows that the band gap of LPPy is lower than 
the one of PPy, as seen in Table [Q and Fig. For 
the LPPy, the dashed line is related to a nearly free 
electron states (NFE). This state is characterized by or- 
bitals located outside the polymer, such states are also 
seen in graphite 30j. The gap for LPPy reported in Ta- 
ble J] neglects the NFE states. Nevertheless, the change 
in the backbone symmetry should generate an increase 
in the band gap even for the LPPy. This suggests that 
another competing phenomenon must be taken into ac- 
count to explain the band gap variation of all ladder-type 
polymers. 

In order to find the difference between the pyrrole and 
thiophene polymer families, the band structures for their 
fully relaxed backbones only was computed. These corre- 
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spond to the cis-polyacetylene for LPT and LPPy while 
for PT and PPy it is a mix between cis and trans- 
polyacetylene. The impact on the band gap of adding 
a sulfur or nitrogen atoms to these backbones can then 
be studied. Within LDA, the gap for cis-polyacetylene is 
0.5 eV, compared to the 1.58 eV and 1.61 eV for LPT 
and LPPy respectively. This band gap increase is sim- 
ilar for the two polymers, which demonstrates that the 
effect of the added atoms is important but also that it 
doesn't depend strongly on its nature. The bond alterna- 
tions in these ladder-type polymers are similar and thus 
their band gaps have roughly the same value. 

For the non-ladder polymers, the backbone presents 
a LDA gap of 0.8 eV, in contrast with 1.38 eV for PT 
and 1.76 eV for PPy. A difference between adding an 
nitrogen or sulfur atom is now observed. The ~0.4 eV 
increase of the band gap between PPy and PT can't be 
accounted for by the variation of the bond length alter- 
nation as seen in Fig. ^ Thus, the intrinsic properties 
of the added atoms must be considered to explain the 
difference in band gaps. The larger electronegativity of 
nitrogen as compared to sulfur causes a greater electron 
attraction in its local environment. This implies that its 
nearest neighbor atoms should be positively charged. At 
a smaller scale, the same phenomenon is found in PT. As 
mentionned before, these polymers are composed of two 
types of carbon atoms. In the case of PPy, the charge 
disparity between the two types of carbon is 0.3e, in con- 
trast with PT which is 0.17e. This results in different 
hopping matrix element parameters in a tight-binding 
model. Since the backbone behave more like a trans- 
polyacetylene, there will be a larger opening of the gap 
for the PPy than for the PT. Hence, to fully explain the 
electronic properties of these polymers it is necessary to 
examine the charges distribution in the polymer's back- 
bone. 

In the case of the ladder polymers, the sulfur or nitro- 
gen atoms still attract electrons in their surroundings, as 
seen in the Bader charge values reported on Fig.Q] But, 
in the LPT and LPPy, all carbon atoms are identical 
by symmetry. Therefore, all the charges are equal. This 
implies that the charge alternation is not a variable in 
our tight-binding model and would not contribute in the 
band gap value. This explain why it is not necessary 
to take into consideration the charge variation in those 
ladder-type polymers to determine their band gaps. 



CONCLUSION 

To conclude, the present study fully explains the be- 
havior of band gaps for two groups of polymers. It was 
shown that the band gap of the LPT was larger com- 
pared to PT. For the pyrrole family, the results found 
were opposite, the band gap of LPPy being lower than 
for PPy. To account for these results, a competition be- 



tween two phenomenons is presented, the bond alterna- 
tion of the atomic structure and the charge of the carbon 
atoms. These two, depending which one dominates, can 
give rise to an increase or decrease of the band gap of the 
polymer in consideration. 
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